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We discuss the influence of lowest Landau level (LLL) fluctuations near H C 2(T) on flux lattice 
melting in YBa2Cu307-s (YBCO). We show that the specific heat step of the flux lattice melting 
transition in YBCO single crystals can be attributed largely to the degrees of freedom associated 
with LLL fluctuations. These degrees of freedom have already been shown to account for most of 
the latent heat. We also show that these results are a consequence of the correspondence between 
flux lattice melting and the onset of LLL fluctuations. 
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Recent high-quality specific heat measurements |l| |j] 
on YBa2Cu307_ d - (YBCO) single crystals have uncov- 
ered features that leave little doubt that the flux lattice 
melts via a first order phase transition into a new state 
of matter called a vortex liquid. B The sharp spikes ob- 
served in these measurements at fields of up to 16 Tesla 
(T) reinforce the previous indications of flux lattice melt- 
ing (FLM) in earlier specific heat measurements [&-[?J and 
resistivity and magnetization experiments. [|8| — [lOfl There 
are many theoretical ]Tl| — |l6|| treatments and numerical 
simulations |l7|-|2^| which treat the flux lattice melting 
in both the low and high field regimes. The considerable 
breadth of the theoretical approaches brought to bear on 
the question has produced a considerable number of in- 
sights, but no clear overall consensus has yet evolved as 
to the origin of the features that the experiments have 
revealed. 

A prominent feature of the specific heat results is the 
spikes associated with the heat of melting. Along with 
these spikes, steps were reported in fields up to 9T in 
Ref. S, with a larger specific heat on the vortex liquid 
side of the transition. Such steps were also observed in 
Refs. 0Jm|. Schilling et al. || found that they were un- 
able to explain the steps in terms of the Abrikosov ratio, 
the effective Debye temperature, the number of vortices 
or the vortex degrees of freedom. Indeed, other authors 
have also shown that degrees of freedom not associated 
with the vortices contribute a significant amount to the 
entropy jump (i.e. the specific heat spike) at the FLM 
transition. For example, Hu and MacDonald found 
in their Monte Carlo study that 90% of the latent heat at 
the FLM transition comes primarily from "the change in 
entropy content at microscopic length scales associated 
with a change in the magnitude of the superconducting 
order parameter and not from changes in the entropy con- 
tent of vortex configurations." |21| Alternatively, it was 
suggested by Volovik that some of the excess entropy 
could be attributed to "electronic" degrees of freedom in 
the vortex background, that is, to quasiparticle excita- 
tions close to the gap nodes of a d-wave superconductor. 



In any case, it appears that the vortices are not the lead- 
ing contributors to the latent heat. 

In this paper, we address the question of the specific 
heat step which is observed in connection with the spike. 
We show that, for fields larger than 2-3 T, the entropy 
from lowest-Landau-level (LLL) fluctuations provides a 
significant, if not leading, contribution to the steps ob- 
served in the specific heat at the FLM transition. Such 
an explanation for the steps implies a deeper connection 
between FLM and LLL fluctuations: namely, that flux- 
lattice melting corresponds with the onset of LLL fluc- 
tuations. Evidence for such a correspondence has been 
presented before, p3] but the arguments for this idea 
will be reinforced and made more persuasive here. We 
develop our argument along the following lines: First, we 
will show that an analysis of the specific data of Ref. B] 
in terms of the analytical LLL expressions of Refs. p5||26[ 
indicates that a good portion of the "step" near the spike 
can be attributed to the onset of LLL fluctuations. Sec- 
ond, a comparison of the spike positions with LLL pre- 
dictions |l5|,^7],^8| emphasizes the correspondence of flux- 
lattice melting with the onset of LLL fluctuations. 

It is helpful to briefly review fluctuations in super- 
conductors and the terminologies commonly employed 
to describe them. We are particularly concerned here 
with LLL fluctuations. In conventional bulk supercon- 
ductors, there is a phase transition at H c2 (T). Whether 
or not fluctuations are important is determined by the 
Ginzburg criterion. [ p9| If fluctuations are negligible, the 
signature of the transition in specific heat measurements 
as a function of temperature consists of ramps with the 
mean-field discontinuity at the transition. pc|| If fluctu- 
ations are significant, they contribute a "bump" on top 
of this mean-field ramp and there is no sharp discontinu- 
ity. |3l| These fluctuations might be generally denoted as 
H C 2 fluctuations. At higher fields (larger than about 1- 
2T as specified below) , the system can be treated within 
a Ginzburg-Landau, lowest-Landau-level formalism. In 
this case, these fluctuations are called LLL fluctuations. 
Fluctuations can contribute to the entropy through mi- 



1 



croscopic order parameter amplitude fluctuations as well 
as vortex position fluctuations, represented by zeroes of 
the order parameter. The former are those found re- 
sponsible in the simulations of Hu and MacDonald pl| ] 
for most of the entropy change. 

Analytical expressions for the specific heat of two- 
dimensional (2D), layered, and three-dimensional (3D) 
superconductors have been derived [^5 26 through the 
use of a Ginzburg-Landau LLL approach. We will use 
these theoretical results to fit the YBCO, H \\ c data 
(for field values 3T-9T) of Ref. §. We will use the 3D 
expressions, |^(| since YBCO is relatively isotropic when 
compared to the Bismuth-, Mercury-, or Thallium- based 
copper oxides. This 3D function is written in detail as 
Eq. (2) of Ref. (24| and it provides a good fit to the data 
for fields H>2T for YBCO type materials. In 
the data we consider from Ref. || , the zero-field data was 
subtracted off. Since 3D LLL expressions do not apply to 
zero field, we will approximate the theoretical zero-field 
contribution by a mean-field expression. This is valid for 
the temperature ranges that we are investigating since 
the fluctuations in zero field are negligible here. 



o 
E 

i- 

<o 
O 

t; 
h- 
x 
o 



15 
5 



-15 



-25 



3T 




4T 




5T 




6T 




71 „iL~— - 


8T 
9T 











74 



78 



82 
T[K] 



86 



90 



FIG. 1. The H || c YBCO specific heat data of Ref. § 
along with fits, discussed in the text, to the 3D specific heat 
function of Ref. §6). 

In Fig. [j], the best fit to Schilling et al.'s data || is 
shown. The fits are quite satisfactory, except of course 
in the regions of the spikes, and the curves provide a 
smooth crossover from the vortex solid phase to the vor- 
tex liquid phase. Thus, we can see that a major portion 
of the "step" near the spike can be attributed to the on- 
set of of H c2 or LLL fluctuations. The remaining part of 
the "steps" might be explained, at least to some extent, 
in terms of the thermodynamic-equilibrium properties of 
the first order vortex lattice phase transition discussed 
in Ref. Q and perhaps also in part by the quasiparticle 
excitations. Q The fitting parameters used in Fig. |l| are 
the c-axis coherence length £ c = 7.218A, the ratio of the 
slope of H c2 (T) to the Ginzburg-Landau parameter k: 



H' c2 /k = 3.40 x 10 2 T/K, the mean-field transition tem- 
peratures T C {H) = 91.33, 90.86, 90.50, 90.25, 89.8, 89.51, 
89. 2K for H = 3T - 9T, and the parameters Q = 8.29, 
K = -1.44, and M = 5.21 of Ref. §§]. The values of 
all parameters are reasonable. For the YBCO materials 
typical values are H' c2 = 1.8T/K, k = 52, and £ c ~ 3A 
The values of T C (H) do produce an H' c2 which is large. 
The most likely reason for this, in our opinion, is the fact 
that the function does not account for 3D/2D crossover. 
An extensive discussion of this point is in Ref. p4| . 
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FIG. 2. The positions of the FLM features of Ref. (| (tri- 
angles) and Ref. [S| (plus signs) plotted in H-T space along 
with their respective fits (see text) to Eq. (|l]). 

We turn now to the second point. Our evidence that 
the LLL fluctuations contribute more to the specific heat 
on the vortex liquid side of the transition than the vortex 
lattice side implies that FLM coincides with the onset of 
LLL fluctuations. To develop this correspondence, we 
turn to the 3D LLL calculation prediction that the melt- 
ing temperature Tm(H) should occur at a fixed value of 
the reduced temperature y = (T U -T C (H))/(T M H) 2 / 3 = 
constant. Her but and Tesanovic [jl5],^7| have calculated 
the value of the scaling constant using density functional 
theory finding, 



(T M - T C (H))/(T M H) 2 ^ 



(327rVi0.5T c0K 2 ^ fc fc B /[<^ 2 £ c 



N2/3 



(1) 



(A similar value for the constant was calculated by 
Hikami et al. p8| using perturbative series expansions.) 
Here we compare the experimental features in the spe- 
cific heat data, which mark the melting, to Eq. (l|). 
We have done this for two sets of data, as shown in 
Fig. ^. The spikes of Ref. || are denoted by the trian- 
gles. The dashed line through them is a two-parameter 
fit of this theory to the positions of the spikes, using 
H' c2 = 1.8T/K, and a linear T C {H). We find a mean-field 
transition temperature T c o = 93.07-fT and that the con- 
stant in the above equation equals to 0.1379-ftT 1 / 3 /T 2 / 3 . 
The standard deviation is 0.05. The value of the con- 
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stant can be calculated from the right hand side of the 
equation using H' c2 = 1.8T/K, k = 52, £ c = 3 A, and 
£ a b = 17.8A These are all within reasonable range. We 
have done a similar analysis for the features observed in 
Ref. (plus signs in the Figure) associated with FLM. 
This fit is also shown in Fig. |. We find T c0 = 92.92^ 
and const=0.1427if l / 3 /T 2 / 3 which would correspond to 
t; a b = 18.26A. As one can see in Fig. |^, the fits to both 
sets of data are very good. There is somewhat more devi- 
ation at the lower fields (H ~ 1 — 2T), which is reasonable 
since that is where the LLL approximation is expected 
to break down. That the values of the fitting parameters 
to data from two separate YBCO samples are reason- 
able and nearly the same reinforces the idea that FLM 
corresponds with the onset of LLL fluctuations. 

Evidence for the correspondence of FLM with the onset 
of LLL fluctuations has been previously found 0] using 
the approach of Roulin et al. [|| These authors identi- 
fied the peaks in the differential C{H + SH, T) - C(H^ T) 
with the flux lattice melting temperature. In Ref. |24| , it 
was shown that the peaks in the differential could be par- 
tially accounted for by the onset of LLL fluctuations. In 
particular, peaks in the differential of "theoretical" data, 
generated using the functions of Ref. ^(|, were used to 
identify the temperatures of the onset of LLL fluctuations 
which were then shown to correspond with flux lattice 
melting temperatures found in experiments. [|||] 

Associating FLM with the onset of LLL fluctuations 
may have escaped previous researchers because one does 
not expect such fluctuations to extend to temperatures so 
much lower than T c2 (-ff). Yet, simple estimates using the 
Ginzburg number do reveal that in zero field fluctuations 
can become significant at temperatures on the order of 
five Kelvin Jl~5| ] (even more in the presence of a magnetic 
field) below this temperature. Furthermore such a corre- 
spondence is not inconsistent with the behavior of con- 
ventional superconductors where FLM and H C 2(T) are 
indistinguishable since the Ginzburg criterion is several 
orders of magnitude smaller than in high-temperature 
superconducting materials. 

The statement that FLM corresponds with the onset 
of LLL fluctuations could be recast in terms of a field- 
dependent Ginzburg number Gi{H). One can simply say 
that FLM is determined by Gi(H). The usual Ginzburg 
number is only defined in zero field. The field-dependent 
Ginzburg criterion says that fluctuations become impor- 
tant when (T - T C (H))/T C (H) ~ Gi(H). Since we have 
found evidence that the fluctuations become important at 
the FLM temperature, Gi(H) could then be introduced 
from Eq. pi). One finds, 



Gi(H) ~ H 2 / 3 (32TT 2 VT0^K 2 ^ ab k B /[ 



^ 2 £c]) 2/3 . 



(2) 



observed at the FLM transition in the high-quality spe- 
cific heat measurements of Schilling et al. |}| originate 
mainly in the entropy associated with lowest Landau 
level fluctuations. Thus, the step appears to be amenable 
to the same explanations as those for the large entropy 
jumps given in Ref. |2l]]. We have further shown that 
the FLM features correspond with the onset of LLL fluc- 
tuations and have derived an expression for the field- 
dependent Ginzburg number that applies to fields where 
the LLL approximation is valid. We speculate that at 
lower fields, FLM corresponds with the onset of what we 
have called H C 2 fluctuations. 
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This value is seven times larger than the estimate given 
by Blatter et al. Q 

In summary, we have shown that the specific heat steps 
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